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The two dimensional abelian Higgs model is revisited. Using the Dirac formalism for constrained systems, 
it has been established [1] that this model, in its gauge invariant physical sector, corresponds to the coupling 
of a pseudoscalar field, namely the electric field, with a real scalar field, as a matter of fact the Higgs field 
which is the radial component of the original complex scalar field in field configuration space. At the 
classical level, the Euler-Lagrange equations of motion lead to a system of coupled non-linear equations of 
which the spectrum of solutions includes solitons. The linearised spectrum of fluctuations of these solitonic 
solutions has been identified in order to ascertain their classical stability. 

This contribution is organised as follows. Section [2] discusses the Euler-Lagrange equations of 
motion. Upon compactification of the spatial dimension into a circle, in Section [3] static solutions to these 
equations are constructed in closed analytic form in terms of the Jacobi elliptic functions. In Section f?] the 
linearised spectrum of fluctuations for these classical solutions is identified. Some concluding remarks are 
provided in Section [5l 



In this paper the two dimensional abelian Higgs model is revisited. We show that in the physical 
sector, the solutions to the Euler-Lagrange equations include solitons. 



1 Introduction 



1 Contribution to the Proceedings of the Fifth International Workshop on Contemporary Problems in Mathematical 
Physics, Cotonou, Republic of Benin, October 27-Novembcr 2, 2007, eds. Jan Govaerts and M. Norbert Hounkonnou 
(International Chair in Mathematical Physics and Applications, ICMPA-UNESCO, Cotonou, Republic of Benin, 2008), 
pp. 164-169. 
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2 The Euler— Lagrange Equations of Motion 



£ P hys = - -e 2 B 2 + -{d^pf - V{p) - do{-BdoB} + & [B(d W + eA )} , (1) 



The 1 + 1 dimensional abelian Higgs model is described by the Lagrangian density 
£ = - \ [d fl Au - d v A„] [d»A v - d v A») + \{dfj, + ieA„) 0| 2 - V(\<f>\), 

where A^ is the gauge field with gauge coupling constant e and <f> a complex scalar field with self-interactions 
described by the U(l) gauge invariant potential V(|0|). Choosing a parametrisation of the complex field 
as 4> = pe ltp /y/2 and through the Dirac formalism for constrained systems [2] the model in its physical 
sector is described by [1] 

\^{d,Bf \e 2 B 2 + l -{d,pf V{p) d [j 2 

where B = — \E and E is the electric field. The Euler-Lagrange equations of motion for p and B are, 
respectively, 

3>+^(^B) 2 + ^ = 0, d p (±d»B^+e 2 B = Q. (2) 

These equations are non-linear and coupled, and not all their solutions can be analytically found. However 
it might be possible to find analytical solutions for static states of finite energy. In order to reduce 
the difficulties of solving these equations, we consider configurations where the electric field vanishes, 
E(t, x) = 0. In fact the total energy is E p i iys = + E p , where E^ is the total kinetic energy and E p the 
total potential energy of the fields, with, 

£ *=/>(^4 B)2+ ^> 2 

and 

From these expressions, one notices that if B is non zero but p vanishes for some value of x, we have a 
singularity in the equations of motion and furthermore the energy may become infinite, unless the electric 
field vanishes faster than p in such a manner that the quotient of these zeros remains finite. But such 
behaviour induces large spatial gradients in both p and B, implying a large value for the total energy. Hence 
in order to minimise the energy, one needs to consider configurations with B = 0. Another advantage of 
this restriction is that the equations are no longer coupled and we simply have to solve 

r)V 

B = 0, % p+ —=0. (3) 

A further assumption is required though, since it is still difficult to construct all the time dependent 
solutions to the above equation (except for those obtained by Lorentz boosts from a static solution). 
Moreover, any time dependence for a solution increases its total energy. Consequently we restrict further 
to static configurations in which p only depends on the space variable x, leading to the single non trivial 
equation, 

5(M) = ' d^T ' (4) 

where for the potential energy henceforth the Higgs choice will be made, V(p) = (1/8)M 2 (p 2 — p 2 ,) 2 , with 
M > a mass scale and po the expectation value of the scalar field. 



3 Static Solutions 

In order to solve the above equation, let us choose to compactify space into a circle of length 2L with 
—L < x < L. For any constant values of p associated to the minima of the potential, where the potential 
also vanishes given our choice of subtraction constant, the non-linear equation in Q is satisfied. These 
constant solutions B{t, x) = 0, p(t, x) = ±po correspond to vacuum configurations, of which the total 
energy values are minimal and vanish. 
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Besides these configurations, we also have non constant static solutions which behave like solitons. 
Given the equation for p, 



P 



dV 
~dp 



0. 



(5) 



where p"{x) stands for the double derivative with respect to x, there exists a conservation law expressed 
as 



1 



((pi(x)) 2 = V( P )+C , 



(6) 



Co being some integration constant. In this form the solution is readily obtained by quadrature. On the 
circle of length 2L one finds, 



Ps{x) = ±p 



2k 2 



1 



1 



-Mp 



1 + fc 2 ' 



it being understood that the solutions obey either periodic or anti-periodic boundary conditions given the 
remaining freedom in the sign of p existing for the choice of polar paramctrisation of the complex scalar 
field (f> [1]. These solutions are thus expressed in terms of the Jacobi elliptic functions [3]. According to 
whether a periodic or an anti-periodic solution is obtained, the remaining single integration constant is 
related to the elliptic modulus k through the condition, 



LMp 



A(n + r)K(k) 



1 



n 6 N, 



: 
1/2 



periodic, 
antipcriodic. 



Their total energy is represented through the expression 

1 



E 



phys 



dx\-M 2 (p 2 s ~p 2 ) +C 



By substituting the above explicit expression, one finds 



E 



phys 



LM*p% 



2(n + r)K{k) J 



2{n+r)K(k) 



dy 



(1 + fc 2 ) 2 



2k 2 



(1 + fc 2 ) 



sn 2 y - 



2fc 4 



(1 + fc 2 )' 



4 

sn y 



This quantity is finite since < fc 2 < 1 et < sn 2 y < 1. This very fact together with the spatial profile 
p(x) of these configurations justifies their interpretation as solitons. In particular the periodic solution 
with n = 1 in the decompactification limit L — > oo reduces to the celebrated kink solutions of the </> 4 real 
scalar field theory in 1+1 dimensions. 



4 Spectrum of Fluctuations 

Given the explicit solutions of the previous Section, in the present Section we address the issue of the 
classical stability under linearised fluctuations in field configuration space. This requires the computation 
of the spectrum of fluctuation eigenvalues, to ascertain that none of these is negative, which otherwise 
would establish that some modes have an unbounded above exponentially growing amplitude, spelling 
disaster for the corresponding solution. 

Let us consider arbitrary time- and space-dependent fluctuations around the identified solutions, 
B(t,x) = 8B(t,x), p(t,x) = p s {x) + 5p(t,x). The corresponding linearised Lagrangian density, expanded 
to second order in these fluctuations, is 

C = ^ (d»5B) (d»6B) - \e 2 (SB) 2 + \ (d^) (d^p s ) + (d„p s ) (d^5p) 
+ i (d M Sp) (d^Sp) - V( Ps ) - SpV'( Ps ) - \ (Sp) 2 V"{ Ps ) . 
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Applying the variational principle in the fields SB and Sp, the linearised Euler-Lagrange equations for 
these fields read 

" 1 



d M d"6p + V"(p a )6p = 0, 



e 2 5B = 0. 



(7) 



These equations being linear in the field variables, a normal mode expansion is warranted, with a spectrum 
of eigenfrequencies to be identified through the following ansatz for time dependence (the general solution 
then being constructed through linear combinations), 



8p(t,x)=f(x)e-^ t + r(x)e i 



6B(x, t) = g{x)e' l ^ t + g*(x)e luJ2t . 



(8) 



If there exist solutions for which either of the eigenfrequencies u>i or ll>2 is pure imaginary, namely such 
that Lof < (i = 1,2), this would imply an exponential run-away time dependence for at least one of 
the linearised fluctuation modes, hence instability of the corresponding classical solution. In the above 
parametrisation of the normal modes, f(x) and g(x) stand for complex-valued functions on the circle, with 
g{x) periodic and f{x) periodic or anti-periodic according to the periodicity properties of the classical 
solution of which the stability needs to be established. 

By direct substitution of the above ansatz into ©, the following eigenvalue equations are derived, 



d 



dx 2 
1 d 



+ V'XpS} f(x) = ojI f(x) , 



dx p 2 (x) dx 



e 2 p 2 s (x) 



g{x) = lo\ g{x) 



(9) 
(10) 



In the case of the vacuum configuration, B = 0, p = ±poi the solutions to ([9]) and (jT0|) are readily 
constructed through a Fourier scries analysis on the circle for the unknown functions g and /, with 
g(x) = ^/C2L)J2n=-oo et2: ^ x 9n an d f{x) = 1/(2L) Y^~t=-oo e% ^ x 9n- The spectrum is then found to be 
given as 



fn 
9n 



Till 

T 



V L 



M 2 p 2 , 



2 2 

e P , 



which is indeed positive definite. As it should, the vacuum configuration is indeed stable against all possible 
fluctuations in the fields. 



For the non-trivial soliton configurations, the eigenvalue problem reads, 



dy 2 



6k 2 sn 2 y 



2uj 2 



(i + fc 2 )/(y), 



(ii) 



tp(y) 



2(1 + k 2 )ul 
M 2 pI 



where 



W B (y) = 2 



2e 2 k 2 , 1 
=— sn ii 

M 2 y 



sn 2 y 



(12) 



(13) 



In order to obtain (fT2|) . the following change of variable was introduced in (fTU| . g(x) = p s (x)ip(x). As 
a matter of fact, pip is a Lame equation, of which the solutions have been discussed and classified in 
Ref. [3], with the results listed in the Table below. 

As may be seen from those results, the eigenfunctions f(y) = sny ■ cny and f(y) = sny • dny both 
have a positive eigenspectrum of lo\ values and are both antiperiodic. However the last three solutions 
in the Table correspond to periodic functions, one of which is the zero mode associated to infinitesimal 
spatial translations, while among the other two there always exists one with a strictly negative eigenvalue 
ui 2 . Consequently the solitons possessing a periodic boundary condition on the circle are unstable against 
fluctuations corresponding to that specific mode. 
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f(y) 




sn y ■ cn y 


1 M 2 „2 3 


sn y ■ dn y 


1 M 2 2 3fc 2 


cn y ■ dn y 









sn 2 y - ^(1 + fc 2 ± - ^ 2 (1 - fc 2 ) 



For what concerns the function g{x) and its fluctuation spectral equation (|12[) . the latter eigenvalue 
problem being equivalent to solving the Schrodingcr equation for the potential Wb{u), the issue may also 
be addressed from that point of view. Even though the Schrodinger equation may not be explicitly solve 
for that choice of potential, using the fact that the potential (|13[) is positive, implies in any case that the 
spectrum of uj\ eigenvalues is likewise positive. No instability may arise in the g(x), namely the B sector 
of fluctuations. 

5 Concluding Remarks 

The vacuum configuration has thus been confirmed to be stable against all fluctuations. Non-trivial static 
soliton configurations are also stable, but only in the sector of anti-periodic solitons and thus an anti- 
periodic boundary condition on the function f(x). Periodic solitons, though, are always unstable. 

Even though this work has some overlap with some previous studies, its originality lies with the 
fact that it considers specifically only the physical sector of the 1+1 dimensional abelian U(l) Higgs 
model without applying any gauge fixing procedure whatsoever. In particular the factorisation between 
the actual physical sector and the decoupled gauge variant sector offers two advantages. First, that no 
artefacts due to gauge fixing are introduced. Second, that potential spurious instabilities lying solely within 
the gauge variant and unphysical sector are avoided from the outset, a feature any other approach having 
been developed so far cannot achieve. All these approaches have until now relied on some gauge fixing 
procedure rather than applying such a gauge invariant and physical factorisation. Hence these approaches 
runs the risk of identifying instabilities which in fact are pure gauge and thus unphysical. Such difficulties 
are avoided from the outset within our approach. 

When the potential V(p) is taken to be the Higgs potential, as we did, the static solutions in p with 
B = are in fact those of Refs. [4,5]. These authors also solve the model of a single complex scalar field 
on the circle, while in Ref. [6], the authors consider the 1 + 1 abelian Higgs Model on the circle, but in 
contradistinction to what we have done, they apply a gauge fixing procedure in terms of the gauge field 
components Aq and A\. We recover the same solutions on the circle. Our work has established that the 
soliton configurations with anti-periodic boundary conditions are stable. In Refs. [5,6], the static soliton 
solutions with periodic boundary conditions are also found to be unstable. 

A possible continuation of the present study would be the computation of the complete spectrum 
of fluctuation eigcnfrcqucncies for the stable and unstable classical soliton configurations, inclusive of 
the electromagnetic sector contributions but retaining only the physical degrees of freedom, in order to 
determine the quantum corrections of the mass spectrum of these solitonic field configurations of the 1+1 
dimensional abelian U(l) Higgs model, whether on the circle or on the real line. 
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